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1. Introduction. There is hardly a more basic stochastic model 
£> \ than a random walk, and for random walks with negative drift, 

L a basic issue of study is the distribution of its global maximum. 

(T) ■ One reason for interest in this quantity is its connection to queueing 

processes. For a GI/G/1 queue, which is stable in the sense that the 

vq \ mean interarrival time exceeds the mean service time, the waiting 

<^ ' time that the ra-th arriving customer needs to wait until service 

c| ■ begins has a limiting distribution as n tends to infinity. This is 

given by the distribution of the global maximum of a random walk 
with negative drift; see Asmussen (1987, §111.7). In an insurance-risk 
setting, the distribution of the global maximum of a random walk 
with negative drift directly appears in computing ruin probabilities 
over an infinite horizon; see, for example, Embrechts, Kliippelberg 



> 

X 
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03 ■ and Mikosch (1997, 51.1) 



Let (X n ) n>1 be a sequence of independent and identically dis- 
tributed random variables having negative mean. The associated 
random walk is defined by Sq = and for any integer n pos- 
itive, S n = X\ + ■ ■ ■ + X n . The distribution of its maximum, 
M = max n ^ S n , can be represented as a compound-geometric dis- 
tribution as follows. We first agree that the minimum of the empty 
set is +oo. Then, let r denote the hitting time for the positive 
half-line 

r = min{ n : Si < 0, . . . , S n -i < , S n > } . 



This hitting time may be infinite, but it is finite with probability 

p = P{ T < oo } = 1 - P{ Si < 0, S 2 < 0, . . . } . 

Recall that the first strict ascending ladder height distribution is 
defined by 

F+(x) =P{S T <x,t <oo}. 

Since the random walk has a negative drift, F+ is a defective 
distribution with defect 1 — p. By the Sparre Andersen identity 
(Feller, 1971, §XII.7) and Abel's lemma (Karlin, 1975, §11.5), 



l-exp(-^lp{5 n >0}); 



P . , 

\ * — ' n 

see also Chung (1974). It follows that M has a compound-geometric 
distribution subordinate to the distribution H = p~ 1 F + and with 
subordinator a geometric distribution with parameter p. More 
explicitly and following Feller (1971, §XII.5), writing H* n the n-fold 
convolution of H, the distribution W of M is 

W = {l-p)Y J P n H* n . (1.1) 

This step has replaced the original question of analyzing the 
distribution of global maximum of a random walk with the more 
elementary question of analyzing that of a compound sum, at the 
price, however, of introducing a derived distribution, namely, the 
ascending ladder height distribution, which requires its own analysis. 
In the case of a heavy-tailed step size distribution as prescribed 
by a subexponentiality assumption, Veraverbeke (1977) supplies an 
answer to this question through use of the distributional form of 
the Wiener-Hopf factorization. His result establishes inheritability 
of the subexponential property of the right Wiener-Hopf factor, i.e. 
that factor having mass concentrated in positive half line, from that 
of the underlying distribution. To state that result, we agree that for 
any possibly defective distribution function G, we write G its tail, 
that is the function whose value at x is 

G(x) = lim G(t) - G(x) . 

t — >oo 

Consider the Wiener-Hopf factorization F = F + + F_ — F + * 
F_, where F_ and F+ are concentrated on (— oo,0] and (0, oo) 



respectively (see Feller, 1971, §XII.3). Let /j, be the mean of F, 
which we assume to be negative. Veraverbeke (1977) shows that, as 
x tends to infinity, 



F + (z)~L_E f°° F(t)dt. 

-V Jx 



We remark that the Wiener-Hopf factors are given by the strict 
ascending ladder height distribution and the weak descending ladder 
height distribution. 

With these two steps in place, a first-order analysis of the distri- 
bution of M may be completed by using a result on tail-area asymp- 
totics for subordinated probability distributions — in this case for 
the subordinator given by a geometric distribution. For example, the 
result in Athreya and Ney (1972, §IV.4) for first-order asymptotics 
of compound subexponential distributions with geometric subordi- 
nator gives the expected result that, under the assumption of F 
subexponential with negative mean fj,, 

1 _ -1 f°° — 

W(x) ~ F+(x) ~ — / F(t) dt , 

1 ~ P V Jx 

as x tends to infinity. This last step of analysing tail areas for subor- 
dinated distributions is a subject of much interest in the literature. 
In the case of subexponential subordinate distributions, we mention 
the paper, Embrechts, Goldie and Veraverbeke (1979), who prove 
that for any integer valued random variable N independent of the 
sequence {Xi) i>x and such that Ez N is analytic at z = 1, 

P{ S N > x } ~ EN P{X l > x} 

as x tends to infinity. 

Finally, we mention that when F is regularly varying with index 
in the range from —1 to —2 and the mean is finite and negative, 
Omey and Willekens (1986) establish a second-order result for the 
tail W; see also Geluk (1992, 1996). 

In a different range of tail heaviness which we will not explore 
in this paper, Feller (1971, §XII.5, Example c) considers the tail 
behavior of W when F has a moment generating function finite 
in a neighborhood of the origin. He shows that if the moment 
generating function of F is 1 at some positive k, and if the number 
P = J xexp(Kx) dF + (x) is finite, then 

W(x) ~ ^Zl e - KX 



as x tends to infinity. If (3 is infinite, the result should be read as 

W = o(e- KX ). 

2. Expansions. In this paper, we are interested in proving higher- 
order tail area asymptotics for the distribution W of the maximum 
M of the random walk. As previously noted, the representation 
of this distribution as a compound-geometric allows a restatement 
of this problem as one of establishing higher-order results for cer- 
tain compound-geometric distributions. We remark that Omey and 
Willekens (1987) provide second-order results for such a compound 
distribution subordinate to a distribution with finite first moment 
which satisfies certain smoothness and regularity conditions, e.g. 
membership in a subclass of subexponential distributions includ- 
ing distributions with regularly varying tails with index of regular 
variation at most — 1. Note however that their results cannot be 
automatically applied to derive second-order behavior for M. This 
is because their result requires certain smoothness conditions which 
would need to be established for the right Wiener-Hopf factor under 
appropriate conditions on the underlying distribution. 

Our approach to this question is to invoke a result from Barbe 
and McCormick (2004) on asymptotic expansions for tail areas of 
compound sums. Their Theorem 4.4.1 provides an m-term expansion 
for tail area, where the number of terms allowed is constrained by 
smoothness and moment conditions on the underlying distribution. 
To that end, we now present a smoothness condition needed for that 
result. 

Definition. A real measurable function f is smoothly varying with 
index —a and order m if it is ultimately m-times continuously 
differentiable and the m-th derivative f^ m ' is regular varying with 
index —a — m. We denote the set of all such functions by SR- at77l . 

All distributions with regularly varying tails used in applications 
are smoothly varying of arbitrary order. Examples include the 
Pareto, Cauchy, Student, Burr and log-gamma distributions. Any 
function smoothly varying in the sense of Bingham, Goldie and 
Teugels (1984, §1.8.1) is smoothly varying of any fixed order. 

The class SR- at7n may be extended to noninteger orders. This 
is useful to present sharp results. To define SR- atU where w is a 
positive real number, we introduce the following notation. For any 



function h, let 

h(t(l-x)) -h(t) 



Ar tA h ) = s[ g n (x)- 



\x\ r h(t) 



Definition. Let to be a positive real number. Write to = m + r 
where m is the integer part of lo and r is in [0,1). A function h is 
smoothly varying of index —a and order lo if it belongs to SR- atm 
and 

limlimsup sup A£ x {h) = . 

<5^0 t^oo o<|a;|<<5 

We write SR- atUJ for the class of all such functions. 

We remark that the spaces SR- ajUJ are nested, for SR- a>r Z> 
SR_ as for r < s. In particular, if lo is positive with integer part 
m and r = lo — m, membership in SR- a ^ is guaranteed by that in 
SR-a^m+i, that is by checking that the m + 1-derivative is regularly 
varying of index —a — m — 1. For further properties of smoothly 
varying functions of finite order, we refer to Barbe and McCormick 
(2004). 

We now introduce the algebraic formalism to express the result. 
To that end let R m [D] denote the ring of real polynomials in 
D modulo the ideal generated by D m+1 . In other words, any 
polynomial in D divisible by D m+1 is set equal to 0. 

For a possibly defective distribution function G with at least k 
moments finite, we write hg,u its k— th moment. Note in particular 
that Hg,o is the total mass of G, equal to 1 if and only if G is not 
defective. 

Definition. The Laplace character of order m of a possibly defective 
distribution G having a finite m-th moment is the element o/R m [D] 
given by 

l-G,m = Y^ ~kl~ fJ ' G ' kBk ' 

0<k<m 

Note that since the map which associates to a measure its fc-th 
moment is linear on its domain, the map G i— ► l~Q )m is linear on its 
domain. 



The backward signed shift S on polymomials in D is defined 
linearly by SD° = and whenever j is a positive integer, SD J = 
-D?" 1 . It maps R m [D] to R m _i[D]. 

We define the inverse of the differentiation on some functions as 
follows. If / is a function regularly varying of index less than —1, 
we set 

/■oo 

D"7(t) = -/ f(x)dx. 



Clearly, DD _1 is the identity on functions which are regularly 
varying of index less than —1, while D _1 D is the identity on the 
smoothly varying functions of negative index and order at least 1. 

We now present our main result. Recall F denotes the step size 
distribution about which we assume its first moment is negative 
and that its right tail is regularly varying of index —a. The strict 
ascending ladder height distribution is F + and W is the distribution 
of M. Let 



sup< r ^ : / \x\ r dF(x) < oo > . 

J — oo J 



This may be less than a if the lower tail of F is heavier than the 
upper one. In the following theorem, it is implicitly supposed that 
n is greater than 1. 

Theorem. Suppose that F is smoothly varying of index —a and 
order 10. Then, for any integer m at least 1 and less than wAaA/t, 
the moments /j>F_,m and fJ-F + ,m-i are finite and 

TF=(l-p)(Id-/. F+ , m _ 1 )- 2 (S/. F _, m )- 1 (D- 1 F) + (Id- m+2 F). 



Remark. As previously mentionned, Laplace characters of order 
TTi — l are elements of the ring R m _i[D]. The inverses (Id — 
L F+jm _i) -2 and (SZ-i?_ )m ) -1 are taken in that ring, multiplied 
together in that ring, and applied to D -1 ^. 

Remark. The result may seem a little mysterious and not so 
explicit at a first glance. However, the computations related to 
Laplace characters can be implemented with a computer algebra 
package. For instance, the following very short Maple code calculates 
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the expansion given in the Theorem. In that code, Fp and Fm stand 
for F + and F_. 

restart; m:=4: mu[Fp,0] :=l-q: 

LFp:=sum(' (-l)~j*mu[Fp, j] *x~j/j ! ' , ' j '=0. .m-1) : 

SLFm:=sum(' (-l)~j*mu[Fm, j + 1] *x~j/(j+l) ! ' ,'j'=0. .m-1) : 

a:=taylor((l-LFp)~(-2) ,x=0,m-l) : 

b:=taylor(SLFm~(-l) ,x=0,m-l) : 

expand (convert (q*taylor(a*b,x=0,m-l) ,polynom)/x) ; 

One simply replaces x and q in the output by D and 1 — p, with the 
convention that 1/x should be replaced by D _1 . For instance, using 
that hf,i = (1 — p)fJ-F-,i, taking m to be 4, we deduce the 3-terms 
expansion 

W = —D^F + -4— ((1 " PW_,2 - ^F + ,lHF-,l)F 

A*F,1 Z I 1 F,1 

+ 12 3 (3(1 ~p) 2 IJ 2 F _ )2 + 12/UF,l(^F_,lA i i ? +,2 - A*F + ,lA*F_,2) 

+ 36//|_ 1 /i| +)1 -2(1 -p)ix Ftl ix F _^\F' 

+ o(Id" 1 F). 

Remark. An important point to mention with regard to the main 
result is that the expansion it provides for the tail distribution W 
is based on the underlying distribution F of the random walk, its 
derivatives and its integrated tail. This is notable, because the 
starting point to obtain this result is that of a tail area expansion 
for a subordinated distribution based on underlying distribution 
given by F + . Since F + is generally unattainable in an explicit 
form, the formulation of our main result is more attractive than 
that which results from a direct application of a tail area result for 
subordinated distributions. It is a comment on the usefulness of 
this algebraic approach that such an improvement is so easily and 
transparently attained compared to the effort to accomplish the same 
goal analytically. We conclude this remark by noting that our proof 
shows that a penultimate expansion based on F + is given by 

W = (l- p)(Id - L F+ , m )- 2 F + + o(ld- m+1 F) (2.1) 

provided m is less than wA(a-l). When m is 1 in the above, we 
obtain a second-order result in agreement with Theorem 2.2 in Omey 
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and Willekens (1987). Comparing this formula with that given in 
our theorem, we see that the latter is slightly less accurate. The 
reason is that, under the assumption of the theorem, replacement 
of F+ with an approximation based on F, its derivatives and its 
integral comes with a one-order lower error bound, viz. 

F + = (SL F _, m )- 1 D- 1 F + o(Id- m+2 F). 



Application. Finally, we present an application to insurance risk. 
To that end, we introduce some notation. Let Ro = x be the initial 
capital of an insurance company. We assume that the claim amounts, 
(A n ) >1 , are independent, with common distribution function L 
having a smoothly varying tail of index —a and order m + 1. We 
also assume that the inter claim times (T n ) n>1 are independent, with 
common distribution K, and independent of the claim amounts. 
Finally, we assume the intensity of the gross risk premium is some 
positive c. The net loss to the company in period n is X n = 
A n — cT n . The sequence (X n ) n>1 is a sequence of independent 
random variables with distribution F(x) = L L(x + ct)dK(t). 
Under the assumptions on L, it follows that F is ultimately m-times 
differentiate and 

F^ m \x) [°°LW(x + ct) ,„^ 

■ ail (t) . 



U m )(x) J L( m )(x) 

This implies that F is smoothly varying of index —a and order to. 
Let ip(Ro) be the probability of eventual ruin given R$. Writing as 
before S n for the random walk with increment Xi and M for its 
maximum, ip(x) = P{M > x}. We follow our established notation 
and set F + and F_ for the strict ascending and weak descending 
ladder height distributions for the random walk S n . We assume that 
X\ has negative expectation. We have following expansion of the 
ruin probability, which obviously follow from the theorem. 

Corollary. Assume that L is smoothly varying of index —a and 
order m + 1. Assume also that fj, = EA\ —cET\ is finite and negative 
and that to less than a. Then, 

^(x) = (l~p)(ld-L F+>m _ 1 )- 2 (SL F _ >m )- 1 B- 1 F(x) 

+ o(x- m+2 F(x)) . 



3. Proof of the theorem. We first show that /Up_ im is indeed 
finite. The distributional form of the Wiener-Hopf factorization 
implies that on the negative half-line, 

F_ = F + F + * F_ . 

Since F + has defect 1 — p, this yields F— ^ F + pF_ , from which we 
deduce F_ ^ (1 — p)~ 1 F. This proves the finiteness of HF-,m- 

We begin the proof of the main part of our theorem by establishing 
a preparatory lemma. 

Lemma 1. Let {Yi)i>\ be a sequence of nonnegative random 
variables, independent and identically distributed with finite and 
positive mean. Let (Z n ) n>0 be their corresponding random walk. 
Furthermore, let f be a regularly varying function of index —(3 less 
than — 1. Then, 

lim — — - V Ef(x + Z n 



Proof. We may assume that / is ultimately positive. Since / is 
regularly varying of negative index, it is asymptotically equivalent 
to a nonincreasing function (Bingham, Goldie and Teugels, §1.5.2). 
Therefore, we can assume without any loss of generality that / is 
nonincreasing. Let 9 be positive and less than the mean of the Y^s. 
We have the trivial bound 

J2 Ef(x + Z n )^J2 f( x + 0n ) + /( x ) zZ P { Z n^0n}. 

n^O n^O n>0 

Since / is regularly varying, 



/>OG 

^f{x + 9n)~ I f(x + 6s)ds. 



as x tends to infinity. The change of variable 9s = xz yields 



X 



f( x + 9s)ds = - f(x(l + z))dz 

(i Jo 

_ xf{x) 
0(13-1)' 
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Let M be such that E(Y 1 A M) is greater than 9. Such M exists 
by monotone convergence of Y\ A M to Y\. Let (Z n vt ) n -^o be the 
random walk associated to the sequence (Yi A M) i>1 . By the Hsu 
and Robbins theorem (see Chow and Teicher, 1988, §10.4), the series 
^2 n>0 P{Z^f ^ n6 } is finite. Since Z^f is at most Z n , this series 
is at least ^2 n>0 P{ Z n ^ n0}, and the latter is finite as well. 
Therefore, since 9 is any positive number less than EY 1 , 

l 7J* v w(x) g Efix +Zn) ^ w^m • 

To obtain a matching lower bound, let now 9 be a number greater 
than 1, and let e be a positive real number. Since / is ultimately 
positive and nonincreasing, ^2 n>0 Ef(x + Z n ) is ultimately at least 

^Ef(x + n{EYi+e))l{\Z n -nEYi\ < en ; 9~ 1 x^nEY 1 ^9x}. 

But if 9~ x x ^ nEY\ ^ #x, as x tends to infinity, 

l + -(£Yi + e)J . 

Moreover, in that range of n, for x large enough, the strong law of 
large numbers implies that P{ \Z n —nEY\\ ^ ne} ^ 1 — e. Therefore, 
S n >o Ef(x + Zn) is ultimately at least 

(l-e)f(x)^2( 1 + -( EY i + e )) I{0 _1 3 < nEY t ^9x} 

n>0 



f 9x/EY 1 p 

(l-e)/(x)/ (l + -(£Y 1 + e )) ds 



Orr/BYi 

(l-e)/(x): 



1 - (3 EY X + e 



s \ -/3+i 

l + -(£Yi+e 

x 



5-1 



x/BYj 



Since # and e are arbitrary, we can make 9 tend to infinity after 
taking the asymptotic equivalent of the lower bound as x tends to 
infinity, proving that 
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as x tends to infinity. ■ 

Note that with unimportant and additional assumptions an alter- 
nate proof of Lemma 1 based on the renewal theorem may be given, 
slightly shorter, but not as direct. To sketch it, write G the distribu- 
tion function of Y^ and consider the renewal function U = X]n>o G* n - 
We see that J2 n >o Ef{x + Z n ) = f f dU. When / is smooth, an in- 
tegration by parts and a change of variable bring this integral to the 
form s Jj 00 f'{xs)U(x{s — 1)) ds. The renewal theorem (Feller, 1971, 
§XI.3) yields U(x(s — 1)) ~ x(s — 1)/EY\ as x tends to infinity, uni- 
formly in s at least 1. The result then follows by standard arguments 
involving regular variation. 

The main argument for proving our theorem is to show that F + 
is smoothly varying of index — a + 1 and same order as F. This is 
stated in the next lemma. 

Lemma 2. The strict ascending ladder height distribution F + is 
smoothly varying of index — a + 1 and same order to as F. Moreover, 

F + ~ ^ IdF. 

(a - l)fJ,F-,i 

Proof. The proof has four steps. 

Step 1. A representation for F + . By the distributional form of 

Wiener-Hopf factorization, we have 

F = F + + F_ -F+*F-. (3.1) 

It is convenient to introduce the following integral operator, 

U F _g(t)= f g(t-u)dF_(u). 

J — oo 

As usual, powers of operators are defined inductively. In particular, 
Up_ is the identity and Up_ = Uf_ °Up~ l for any integer n positive. 
On (0, oo), we can write (3.1) as F + = F + F + * F_, which leads to 

F + = F + U F _F + . (3.2) 

By recursion this yields 

F + = J2 Up_F + U F ' +1 F + . 
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Note that F_ cannot be the distribution degenerate at since 
F is assumed to have a negative mean. Let (5 / i) i>1 be a se- 
quence of independent random variables, all with the same distri- 
bution F-, and let (Z n ) n>0 be their random walk (note that the 
signs are changed compared to the previous lemma). Observe that 
So<i<n Up F is nondecreasing in n and that, by dominated conver- 
gence, U F ^ 1 F + {x) = EF + (x — Z n+ i) tends to as n goes to infinity. 
Consequently, we obtain the representation 

F + = Y^U t F F. (3.3) 

Note that combined with Lemma 1, this representation yields Ver- 
averbeke's (1977) theorem asserting that 



F + ~- IdF. (3.4) 

(a- 1)Hf-,i 

Step 2. A representation for F + . Let k be a positive integer at 
most w A (a - 1). Using the mean value theorem, there exists a 
sequence of real numbers, (#n) n >o> nonnegative and at most 1, such 
that 

J2\-( U Fy (k ~ 1] (x + e) - L^_F (fc - 1} (x)) - U^_F {k \x) 

n>0 



Y,\UF^F {k) {x + e n e)-U^_F {k \x) 



sC J2E\F^(x + 9 n e-Z n )-F^ k \x-Z n )\. 

Since the absolute value of a difference is at most the sum of the 
absolute values, Lemma 1 shows that the above series is bounded as 
a function of x and uniformly in e in some interval (0, rf). Moreover, 
every summand tends to as e tends to 0. Therefore, the series 
tends to as e tends to infinity. This proves that 



Tf = Y / U F _F {k) (3.5) 



i>0 



on some neighborhood of infinity. 
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— (k) 

Step 3. F + ' is regularly varying. The asymptotic equivalence in 
(3.4) implies that F + is regularly varying with index —a + 1. Recall 
that k is at most a; A (a— 1). By assumption F^ k ' is regularly varying. 
By representation (3.5) and Lemma 1, F^_' is regularly varying of 
index —a — k + 1. Taking k to be m, that is [lo\, this proves that 
F + is smoothly varying of index —a + 1 and order m. 

Step 4- Concluding the proof of the lemma. Following Barbe and 
McCormick (2004), for a function h define 

K,s( h ) = SU P SU P \ A t,x h \ ■ 
t^T 0<\x\^S 

This quantity is nonincreasing in r and nondecreasing in 5. Using 
representation (3.5), we see that 

F^itil-x^-F^it) 

= Y, E(F^ (t(l -x)- Z n ) - F^(t - Z n j) . 

Consider x in the range [—5, 5]\{ }. Factoring t~Z n in t(l—x) — Z n , 
the n-th summand in the series above is at most 

^-M)> (m) (* " Z n )\ X-z n , tm -z n) F {m) ■ 
Consequently, for \x\ positive and at most 8, 



I l,X -f- I ^ / J "T^\ lrl )l'j.\ 



to 1 ^\t) 



A r f ,F^ m K 



H,8 



It follows from step 3, F™ x Idi^™). It then follows from Lemma 
1 and our assumption on F that 

lim lim sup lAJ^F^I =0, 

5^ t^oo 0< | a ,| <(5 

proving the smooth variation of order uj of F + . m 

Finally, we present a technical lemma of some independent inter- 
est, particularly in the light of Marie's (2000) work. It is needed 
for the proof of our main result. We remark that the result is not 
proved under optimal conditions. 
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Lemma 3. Let (ai)o^i<m be a sequence of real numbers with a® 
different from 0. For any nonnegative integer k at most m, define 
the differential operators P k iB) = ^2 0<i<k a^D*. Let -0 be a function. 
Let f and g be two functions smoothly varying with index —a and 
order at least m satisfying the differential equations 

P m . k (B)B k f = B k g + o(ij>) , k = 0, 1, . . . , m. 

Then, viewing P m (B) tnI ra [D], 

f = P m (B)- i g + o^). 

The lemma may be interpreted as saying that if the functions B k g 
have a generalized asymptotic expansion in the the asymptotic 
scale B k f\ then / has a generalized asymptotic expansion in the 
asymptotic scale B k g. 

Proof. Write b k the fc-th coefficient of P m (D) _1 . Then 

P m {B)- 1 g= Y, b ^ k 9- ( 3 - 6 ) 

In this sum, by assumption, we can replace B k g by P m - k (B)B k f + 
o(V>)- Since P m _ fc (D)D fc = P m (D)D fe inR m [D], the definition of the 
b k and (3.6) yield -P m (D) _1 g = / + o(ip), which is the result. ■ 

We now conclude the proof of the main theorem. Using (3.2) and 
applying Lemma 2 and a variant of Theorem 2.3.1 in Barbe and 
McCormick (2004), we obtain for any nonnegative k at most m, 

p( k ) = pW _ LF _ tm _ k Ff + (i d — F+) 
= SL F _, m _ k B k+1 F + + o(ld- m F + ) . 

By Veraverbeke's (1977) theorem or (3.4), this implies 

SL F _, m _ k B k BF + = B k F + o(ld- m+1 F) . 

Applying Lemma 3, we obtain 

DF+ = (SL^^-'F + o(Id- m+1 F) . 
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Hence, integrating, 

F + = (SL F _ :m y 1 B- 1 F + o(ld- m+2 F) . (3.7) 

Representation (1.1) and Theorem 4.4.1 in Barbe and McCormick 
(2004), upon noting that if N is a random variable with geometric 
distribution with parameter p, then ENL H ~ m is p(l — p)(ld — 
pLH,m)~ 2 in M m [D], yield formula (2.1), that is 

W = (1 - p)(Id - L F+ , m )~ 2 F + + o(ld- m F + ) . 

To obtain the statement of the Theorem, we again use representation 
(1.1) and apply Theorem 4.4.1 in Barbe and McCormick (2004) to 
obtain that if m is less than a Acu, 

W = (l- p)(Id - L F+tm _ 1 )- 2 F + + (Id" m+1 F + ) . 

Then, we use Lemma 2 and (3.7) to conclude. ■ 
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